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Abstract 

An analysis of the neutron EDM and of the electron EDM in minimal N=l 
supergravity unification with two CP violating phases is given. For the neutron the 
analysis includes the complete one loop gluino, chargino, and neutralino exchange 
diagrams for the electric dipole and the chromoelectric dipole operators, and also 
the contribution of the purely gluonic dimension six operator. It is shown that 
there exist significant regions in the six dimensional parameter space of the model 
where cancellations between the gluino and the chargino exchanges reduce the 
electric and the chromoelectric contributions, and further cancellations among 
the electric, the chromoelectric, and the purely gluonic parts lead to a dramatic 
lowering of the neutron EDM sometimes below the electron EDM value. This 
phenomenon gives a new mechanism, i.e., that of internal cancellations, for the 
suppression of the neutron EDM in supersymmetric theories. The cancellation 
mechanism can significantly reduce the severe fine tuning problem associated with 
CP violating phases in SUSY and SUGRA unified models. 



1 Introduction 



Supersymmetric models with softly broken supersymmetry introduce new sources 
of CP violation which contribute to the neutron and the electron electric dipole 
moment (EDM). It is known that the exchange of SUSY particles close to their 
current experimental lower limits and CP violating phases of normal size, i.e. 0(1), 
will lead to the neutron EDM already in excess of the current experimental bound 
of 1.1 X lO^^^ecm fTJ. Two possibilities to resolve this problem have been commonly 
discussed in the literature. The first is that the phases are not 0(1) but rather 
much smaller, i.e., O(10~^ — 10^^) |, @|. However, a small phase constitutes 
a fine tuning unless it arises naturally, e.g. as a loop correction. The second 
possibility is that the phases are 0(1), but the supersymmetric spectrum which 
contributes to the EDMs is heavy i.e. in the several TeV region and perhaps 

out of reach of even the LHC. In this paper we discuss a third possibility, i.e., that 
of internal cancellations among the different components of the neutron EDM. 
We shall show that such cancellations can dramatically reduce the neutron EDM 
without either excessive finetuning of the phases or pushing the SUSY spectrum in 
several TeV mass range. One then finds that the neutron and the electron EDMs 
can satisfy the current experimental bounds |]l], |^ with phases not unduly small 
and a SUSY spectrum which is not unduly heavy at least for low values of tan/5. 

Although there are many analyses of the EDMs in supersjTiimetric theories 
most of these §|, |§, are without radiative breaking of the electro-weak 

symmetry and sometimes neglecting the chargino contribution to the neutron 
EDM P, 0. For the neutron EDM there are two operators other than the electric 
dipole moment operator, which can contribute to the neutron EDM. One of these 
is the color dipole operator and the other is the dimension six purely gluonic op- 
erator considered by Weinberg^. With the exception of the work of ref.[|l^ most 
of the previous analyses 0, H, |^, 0, |lT| do not take into account the contribution 
of the color and of the purely gluonic operators with the presumption that their 
relative contributions to the neutron EDM is small. However, it was pointed out in 
ref . [jl2[ that the contributions of the color and of the purely gluonic operators can 
be comparable to the contribution of the electric dipole operator in a significant 
region of the minimal supergravity parameter space [|1^]. Currently there is some 
confusion in the literature regarding the sign of the gluino exchange term[^, |^. 
Further, there is no analysis aside from that of ref. ||^ which gives the complete 
one loop contribution including the gluino, the chargino and the neutralino ex- 
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changes against which the signs of the relative contributions of the various terms 
can be checked. Because of the sensitive issue of the cancellation between the 
gluino and the chargino exchange diagrams crucial to the analysis of this paper, 
we have redone the full one loop analysis of the EDM with the gluino, the chargino 
and the neutralino exchanges. We compare our results to those of Ref. 0] in 
Appendix A. In our analysis we have made the standard assumption of ignoring 
all the generational mixing of quarks and of squarks. 

The outline of the paper is as follows: In Sec. 2 we give the general fea- 
tures of the minimal supergravity unified model and discuss the new CP violating 
phases it supports. In Sec. 3 we give our evaluation of the gluino, the chargino and 
the neutralino contributions to the electric dipole operator. In Sec. 4 we display 
our evaluation of the chromoelectric and of the purely gluonic operator contribu- 
tions. Numerical analysis of results and the phenomenon of cancellation among 
the various components is discussed in Sec. 5. Conclusions are given in Sec. 6. Di- 
agonalization of the squark and of the chargino mass matrices paying attention to 
the phases is given in Appendices A and B along with a comparison of our results 
with those of the previous analyses. 

2 N=l Supergravity and CP Violating Phases 

The analysis of this paper is based on = 1 supergravity grand unified theory 
in which supersymmetry is broken spontaneously via gravitational interactions in 



the hidden sector |jT3|, [1^, and the electro- weak symmetry is broken via radiative 
effects. We assume that the grand unified theory (GUT) group G breaks at the 
scale Mg to the Standard Model gauge group, and after breaking of supersymmetry 
the tree effective theory can be characterized by the following symmetry breaking 
sector at the GUT scale 

Vsb{^) = mlzaz'' + (moAoW^(=^) + B^W^^^ + H.c.) + ]-m.W (1) 

Here W'^'^^ = fioHiH2, with Hi, H2 being the two Higgs doublets, W'^^^ is the su- 
perpotential cubic in the fields, mo is the universal scalar mass, Aq is the universal 

trilinear coupling, Bq is the universal bilinear coupling, and mi is the universal 

2 

gaugino mass. In general Aq, Bq^q, and mi are complex. However, not all the 
phases are physical. It is possible to remove the phases of mo,mi and make BqHq 
real by redefinition of the fields and by doing R-transformations on them. We are 
then left with only two independent phases at the GUT scale. One may choose 
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one of these to be the phase of /xq (^mo), and the other to be the phase of Aq (oao)- 
Using renormahzation group evolution one can evolve ^b(O) to low energy and 
one finds 

VsB = ml\H^\^ + ml\H2\^ - [B^ieijH{Hi + H.c] 

+M^[uluL + dldi] + Mluy[LR + Mld*Jii 
+Ml[9:i), + eleL] + Mle*^en 

v2miy cosp cosp smp 

+^ + maW^ V" + miBB] + AVsb (2) 

where {II, Ql) are the SU(2) (slepton, squark) doublets (the generation indices are 
suppressed), and AVsb is the one loop contribution to the effective potential ITS |. 



In our analysis the electroweak symmetry is broken by radiative effects which 
allows one to determine the magnitude of /iq by fixing Mz and to find the magni- 
tude of Bq in terms of tan P=^^j- In the analysis we use one- loop renormalization 
group equations (RGEs) for the evolution of the soft SUSY breaking parameters 
and for the parameter fi, and two-loop RGEs for the gauge and Yukawa couplings. 
The equations for the gauge and the Yukawa couplings, and the diagonal elements 
of the sfermion masses and gaugino masses are such that they are entirely real, 
while the phase of fi doesn't run because it cancels out of the one loop renormal- 
ization group equation of /x. However, both the magnitudes and the phases of Ai 
do evolve and one has 

dAi- , 16 _ ms _ mo 13 mi , , , . ^ 

-7^ = -(^«3— + 3a2— + T7«i— + 6F*A + (3) 
dt 3 mo mo 15 mo 



dt 3 mo mo 15 mo 



dA^ 
~dt 
dAu 
dt 

dAd 
dt 



_ m2 9 _ mi 
—( 

rriQ 5 mo 



.16 ^ mg m2 13. mi t . 

—as h 3^2 h —ai h 3Y At), 

3 mo mo 15 mo 



,16. mg m2 7 . mi f, 

.^"3 ^ 3^2 h —ai h 3Y Ab), 

3 mo mo 15 mo 



(5) 
(6) 

(7) 
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and 



= - 3(52— + ^ai— + , 8 

at mo 5 mo 



where ai=j^, Y(^y,,d,e)= '^''(An)^ ' where g', are the gauge couphngs and 
are the Yukawa couphngs, and t — In-^. The supergravity model with CP viola- 
tion is then completely parametrized by just six quantities: mo, mi , ^o; tan /3, O^q 
and a AO- There are 32 new particles in this model. Their masses and interactions 
are determined by the six parameters above. Thus the model is very predictive. 

3 EDM Calculation 

One defines the EDM of a spin-| particle by the effective lagrangian 

= —dfi^a^.^^il^F'^'' (9) 

which in the non-relativistic limit gives £/ = dfip^a.EipA where ip^ is the large 
component of Dirac field. In renormalizable theories the effective lagrangian (9) 
is induced at the loop level if the theory contains a source of CP violation at the 
tree level. For a theory of fermion ipf interacting with other heavy fermions ipi^s 
and heavy scalars (p^s with masses m^, m^ and charges Qi, Qk the interaction that 
contains CP violation in general is given by 

- Ant = XlV'/(-f^ife— ^ + Lik^—^)A<f>k + H.c. (10) 

Here JC violates CP invariance iff lm{KikL*j^) ^ 0. The one loop EDM of the 
fermion f in this case is given by 

where A(r) and B(r) are defined by 

1 , 2/nr 
2(i37p^3-r+- 



-4(0 = .,, ^,2 (3-^+^) (12) 



and 

, 1 , 2rlnr, , , 

where one has charge conserved at the vertices, i.e., Qk — Qf — Qi- The loop 
diagrams corresponding to the term A is Fig. (la) and to the term B is Fig. (lb). 
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3.1 Gluino Contribution 

The quark-squark-gluino interaction is given by |T1 



- = V2g,T^, i-Qi^-Y^9ain + qi^A^Ml) + H.c, (14) 

i=u,d 

where a = 1 — 8 are the gluino color indices, and j, A; = 1 — 3 are the quark and 
squark color indices. The scalar fields qi and are in general linear combina- 
tions of the mass eigenstates which are given by diagonalizing the squark {massY 
matrices for u and d at the electroweak scale [ll4| : 



Ml 



' Ml + m„2 + (i - g„ sin^ Ow) cos 2(3 m„(A;mo - /i cot (3) 

(15) 



" V m„(v4„mo - /i*cot/3) M? + m^^ _^ Q„ sin^ cos 2/3 



and 



M? + m^^ - M2(i + Qdsin2 6'vy) cos 2/5 md{A}mo - /itan/3) 



V TTidiAdmo - fi* tan (3) Ml + + M'^Q a sin^ 9 w cos 2/3 



(16) 

where Qm = | and Qd = — |. We note that the and Ad are not independent but 
evolve from the same common Aq at the GUT scale. Further, 0^(the phase of fi), 
aA„(the phase of Au), and (the phase of Ad) are related to just the two phases 
6^^ and at the GUT scale by renormalization group evolution. We diagonalize 
the squark matrices so that 

qL = Dqiiqi + Dqi2q2 (17) 

qR = Dq2iqi + Dq22q2- (18) 

where qi and q2 are the mass eigenstates. A more detailed discussion of the diag- 
onalization is given in Appendix A. 

In terms of the mass eigenstates qi and q2 the gluino contribution is given by 

d!-,iuUe = E Im(rf )5q,B(^), (19) 

where F^^ = Dq2kD*if., as=|^, irig is the gluino mass, and e is the positron charge. 
3.2 Neutralino Contribution 

In order to discuss the neutralino exchange contributions we first exhibit the neu- 
tralino mass matrix: 
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/ fhi —M^ sin 6^ cos P sin 6w sin [3 \ 

^ _ 7712 cos 6*1^ cos /5 — cos sin [3 

^ —Mj singly cos/? cos cos /3 — /i 

\ sin 6*1^ sin /3 — cos6'vi/ sin/? — /i / 

(20) 

The matrix M^o is a complex non hermitian and symmetric matrix, which can be 
diagonahzed using a unitary matrix X such that 

X'^ MyAiX = diag(m^o, m^o, m^o, m^o). (21) 

By rearranging the fermion-sfermion-neutrahno interaction |T2|, the neutrahno ex- 
change contribution to the fermion EDM is given by 



^-n^^ole = E E I-(^/-)#Q/B(^), (22) 

^/i bm C7VI/ ivij:^ ivij^ 

where 

r]fik = [-V2{tan9w{Qf - T-if)Xii + T^fX2i]D*f^f. + KfXuD^j^k] 

{V2t&n9wQfXiiDf2k - KfXuDfik). (23) 

Here we have 

_ rriu _ md,e . .^ 

V 2mvi/ sm /) v 2mvi/ cos /? 



where 6 = 3(4) for T^f = 



-H- 



3.3 Chargino Contribution 

To discuss the contribution of the chargino exchanges we exhibit first the chargino 
mass matrix 

f rh, v/2m^sin/?\ 
Vv2mpycos/) fj, J 

This matrix can be diagonahzed by the biunitary transformation 

U*McV-^ = diag(m^+, m^+) (26) 

where U and V are unitary matrices(see Appendix B). By looking at the fermion- 
sfermion-chargino interaction we find the chargino contribution to the EDMs for 
the up quark, the down quark and for the electron as follows 

e.w,.no/e = E Eini(r„.,)-^[g,B(-^) + (q„ - q,)a(^)]. 



y k=l 1=1 dk dk dk 



(27) 
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Ccw,.no/e = . E E ini(r...)-^[g.B(-^) + (q. - g.)A(-^)], 

47r sm Uw k=i i=i ^^kk J-^kk ^^-^uk 

(28) 

and _ ^ 

dl.Har,^no/e = ^ E m^+Im([/- V^t) A(^) (29) 

where 

Tmfc = KuV*2Ddlk{U*iD*i_]_,^ - KdU*2D^2k) (30) 



Tdifc = KdU*2DuikiVilDl^k - i^uV*2Dl2k)- (31) 

The sum total of the gluino, the neutrahno and the chargino contributions to the 
EDM gives us the total EDM. To obtain the neutron EDM contribution from the 
electric dipole moment operator we use the non-relativistic SU{Q) quark model 
which gives 

dn = ^[Add-du]. (32) 

The analysis of dn above is at the electro-weak scale and it must be evolved 
down to the hadronic scale via renormalization group evolution to give 

d^ = v'^dn (33) 

where rj^ is the QCD correction factor and we estimate it to be 1.53 in agreement 



with the analysis of ref. |jTO| 



4 The Chromoelectric and the CP Violating Purely 
Gluonic Dimension Six Operators 

The quark chromoelectric dipole moment is defined to be the factor d^ in the 
effective operator 

Cj = -'-d'^qa.^^.T-qG^^^ (34) 

where are the generators of SU{'i). The gluonic dipole moment d^ is defined 
to be the factor in the effective operator 

C-I = ~ -^d^ f afi^G a^pCp^G ^Xat^^'^" (35) 
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where 



where Ga/ny is the gluon field strength tensor, fa/s-y are the Gell-Mann coefficients, 
and e'^"^'^ is the totally antisymmetric tensor with e^^^^ = +1. An analysis of 
these operators in minimal supergravity with two CP violating phases was given 
in ref.|]T^. We quote the results from that work here. For the chromoelectric 
dipole moment one has three contributions; from the gluino exchange, from the 
neutralino exchange, and from the chargino exchange. These are given by|12|j 

d^-,u.no = ^ E Im(r-)^C(^), (36) 

fe=l Qk Qk 

, 1 , 2rlnr 18/nr, , , 

o(r — 1) 1 — r 1 — r 

dq-neutralino = E z2 ^"^(^9^^) T7#B( — ^) , (38) 

^'^^ k=li=l ^^^qk ^^^qk 

dq-chargino = Hk) J^t^iTlt) ' (^9) 

k=l i=l ^^^qk ^^^qk 

The contribution to the EDM of the quarks can be computed using the naive 
dimensional analysis |jl6| which gives 

7C _ e yc^c 



and 



< = (40) 

where is the QCD correction factor for the color dipole operator. 
For the CP violating dimension six operator 



= -3asmA'lm{Tl')^^^ll{z^,z„Zt) (41) 



a 

where 



'^'-^\z, = {^f (42) 



m-g m-g 



The contribution to dn from dP can be estimated by the naive dimensional analysis ||16|| 
which gives 

= "-^d^^V^ (43) 

where M is the chiral symmetry breaking scale and has the numerical value 1.19 
GeV, and r]'^ is the renormalization group evolution factor of the dimension six 
operator from the electroweak scale down to the hadronic scale. We estimate that 



rf" rf^ r-u 3.4 in agreement with the analysis of ref. |T0[. To get the contributions 
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of the chromoelectric and dimension six operators to EDM we used the reduced 
couphng constant and the naive dimensional analysis. There is another way of 
estimating this contribution for the chromoelectric operator and that is using QCD 
sum rules The use of QCD sum rules rather than of the naive dimensional 
method would not change the conclusions of this paper. 



5 EDM Analysis 

As already stated while the EDM of the neutron has been analysed in many works, 
most of the previous analyses have been carried out within MSSM. Our analysis 
here is in the framework of N=l supergravity and we use radiative breaking of the 
electro- weak symmetry including one loop effective potential terms []15|] to anal- 
yse the EDMs in the six dimensional parameter space of the theory given by 
rriQ, mi , Aq, tan (3,9^0 and a^o- The constraints imposed on the radiative electro- 
weak symmetry breaking include imposition of color and charge conservation, ex- 
perimental lower limit constraints on the sparticle masses from LEP, CDF and 



DO, and the experimental constraints on 6 ^ s + 7 from CLE0||T8[. (Details of 
the analysis are similar to those of Ref. |p.9||.). As mentioned in Sec.l in most of 
the previous analyses in the literature the effects of the chromo-electric and of 
the purely gluonic operators have been assumed small and ignored. As shown in 
ref.|]12[ this is an erroneous assumption as the relative contributions of the elec- 



tric, of the chromo-electric , and of the purely gluonic operators are highly model 
dependent and their ratios can sharply change as one moves in the six dimensional 



parameter space of the model. In fact, it was shown in ref.[0 that contrary to 
the assumptions generally made the contributions of the chromoelectric and of the 
purely gluonic operators can be comparable to and may even exceed the contribu- 
tion of the electric dipole term. Because of the significant contribution that the 
chromo-electric and the purely gluonic operators can make to the neutron EDM 
we include in our analysis all the three contributions, i.e., the electric, the chromo- 
electric and the purely gluonic operator contributions. However, we do not include 
in the analysis the effects induced by the phase in the Kobayashi-Maskawa(KM) 
mass matrix in the renormalization group evolution of the SUSY phases, since 
these induced effects are known to be very small[^ ^ 

One of the important phenomenon we find in our analysis is the possibility of 
destructive interference between the gluino and the chargino exchange diagrams 
for the electric dipole and for the chromoelectric terms. This generally happens 
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when the signs of phases of 6fj,g and are opposite. In this case there is also a 
destructive interference between the fi and the At terms in the purely gluonic part. 
In addition to the above one also finds a further cancellation among the electric, 
the chromoelectric and the purely gluonic parts. Constraints on the theoretical 
analyses are provided by the experimental upper limits on the EDMs. For the 
neutron the current experimental limit is 

d„ < 1.1 X lO^^^ecm (44) 

and for the electron the limit is [§] 

4 < 4.3 X IQ-^^ecm (45) 



For the muon the current experimental upper limit is c?^ < 1.1 x 10~^^ ecm||23||(at 
95% CL). This limit may improve by up to four orders of magnitude in a new 
proposed experiment at the Brookhaven National Laboratory||24]. However, the 



constraints on the supergravity parameter space from the current limits on the 
neutron EDM and on the electron EDM are already much stronger than what 
might emerge from the improved muon EDM experiment. For this reason we 
focus in our analysis on the constraints coming from the neutron EDM and from 
the electron EDM. However, we shall sometimes also display the muon EDM for 
comparison along with the neutron and the electron EDM. 

We begin our discussion with a comparison between the electron and neutron 
EDMs constraints on the two basic parameters of the theory, i.e., mo and mi/2- 
As may be seen from Fig. (2a) the electron EDM falls off with increasing mo. This 
behavior is easily understood from Eqs.(22) and (29) since as tuq increases A{r)/Mj 
and B{r)/Mj decrease. Using the experimental upper limit of Eq.(45) in Fig. (2a) 
one finds for \Ao\ = 1.0, tan/? = 3.0 and aAo = 6*^0 = Jq the following constraints 
on mo: mo > 1320 GeV for mi/2 = 700 GeV, mo > 1420 GeV for mi/2 = 600 GeV, 
and mo > 1520 GeV for mi/2 = 500 GeV. A similar analysis holds for Fig. (2b). 
Here using the experimental upper limit on the neutron EDM of Eq.(44) and for 
the same above parameters one finds: mo > 2500 GeV for mi/2 = 700 GeV, 
mo > 2680 GeV for mi/2 = 600 GeV and mo > 2840 GeV for mi/2 = 500 GeV. 
Thus in this region of the parameter space the upper limit on the neutron EDM 
rather than the upper limit on the electron EDM is the more severe constraint on 
mo. The dependence of de and dn on mi/2 is displayed in Figs. (2c) and (2d). The 
broad maxima for small mi in these graphs arise from an interplay between the 
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factors 171 + ,771^0 and rrin which increase as mi increases, and the functions A(r), 
B{r) and C(r) which decrease as mi increases. By carrying out the same analysis 
as for the mo dependence, one finds here also that the experimental upper limit 
constraint for the neutron EDM is a more severe constraint than the one for the 
electron EDM. 

The dependence of the EDMs on in displayed in Fig. (2c) and on tan/3 
in Fig.(2f). The conventional fine tuning problem can be understood from the 
analysis of Fig. (2c) where the phase 6fj_o must lie in a very small corridor around 
the origin to satisfy the current experimental constraints on the neutron EDM. 
Fig.(2f) shows that the EDMs arc an increasing function of tan/9. This behavior 
can be understood easily for the electron and for the muon EDM since these involve 
a factor of l/cosf3 which increases as tan/9 increases. For the neutron EDM case, 
there are contributions from both the up quark and from the down quark with 
different tan/3 dependences. However, the down quark contribution dominates 
and as Fig.(2f) shows the neutron EDM is still an increasing function of tan/3. 

In the analysis thus far we did not take advantage of the two independent 
phases. To give a comparison of the results arising in the two cases we consider first 
the case of Fig. (3a) where the signs of and 9^^ are both positive. Here we find 
that there are no large internal cancellations within the various components , d^, 
d^, and d^ so these functions do not show any rapidly varying behavior. However, 
d^ and d^ in this case are negative while is positive over the entire 1^40 1 region 
and there is a cancellation among them. In the region of I^qI < 2.5 the cancellation 
is rather small because d^ is relatively small. However, the cancellation becomes 
more significant for 1^40 1 > 2.5 leading to a dip in the total d„ in this region as 
seen in Fig. (3a). 

We consider next the case in Fig. (3b) when the sign of a^o is switched. Here 
each of the individual components d^, d^, and d^ shows a destructive interference 
giving rise to sharp minima as a function of |Ao|. These minima can be understood 
as follows: For the case of d^ and d^ the minima arise as a consequence of destruc- 
tive interference between the gluino exchange and the chargino exchange in the 
one loop diagrams. This illustrates what we have said previously that the chargino 
exchange contributions are as important as the gluino exchange contributions and 
should be included in the analysis contrary to what is often done in the literature. 
The minimum in d^ in Fig. (3b) has a different origin. It can be understood by 
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examining the expression 

Im(rt^) = . {mo\At\ sin at + sin6l^cot/3), (46) 

where 9^ and are the values of and of at the electro-weak scale. Prom 
Eq.(46) we see that the magnitude of Im(rP) depends on the relative sign and the 
relative magnitudes of 9fj_ and of at- Thus a cancellation occurs between the At 
and the /i terms when and at have opposite signs leading to a sharp minimum in 

as a function of \Ao\. Now each of the three terms , d^, and d^, switch sign 
as they pass their zero values. Thus d^ and d^ are negative below their respective 
minima and become positive after crossing them, while d^ is positive below the 
minimum and becomes negative after crossing it. This complex structure now gives 
rise to two distinct minima in the albegraic sum of the three terms, i.e., in the 
total dn as may be seen in Fig. (3b). We pause here to note that for the case when 
there is destructive interference between the gluino and the chargino case, and a 
further cancellation among the electric, the chromoelectric and the purely gluonic 
terms as is the case for Fig. (3b), one finds a drastic reduction in the magnitude of 
dn often by a factor 0(10 - 10=^). 

In Fig. (3c) we give a plot of the EDMs of the electron, the muon and the 
neutron as a function of mi/2 showing the cases when ccaq is positive and when 
aAo is negative. Here for the case when is positive the neutron EDM is large 
enough that it violates the current experimental bound in the entire range of 
"^1/2 < 750 GeV. However, for the case when aAo i^ negative the neutron EDM 
lies below the experimental upper limit for mi/2 > 300 GeV. The large disparity 
between the magnitudes of the neutron EDM for the aAo positive case vs for the 
aAo negative case can shift the balance between which of the two experimental 
constraints, i.e., the experimental upper limit constraint on the neutron EDM or 
the experimental upper limit constraint on the electron EDM, is the more stringent 
one. It can be seen that for the case of constructive interference the experimental 
upper limit constraint on the neutron EDM is generally the more stringent one 
while for the case of destructive interference involving large cancellation it is the 
experimental constraint on the electron EDM which may be the more stringent 
constraint. We shall exhibit this effect further in the analysis of Fig.(3e). 

The destructive interference between the different contributions exhibited in 
Figs. (3a) - (3c) is not an isolated phenomenon but rather a common occurrence 
in a large part of the parameter space. Thus, cancellations occur naturally over 
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the entire parameter space with the appropriate choice for the relative sign of 
and a AO- Further, these cancellations can become exceptionally large in certain 
regions of the parameter space. An example of this effect already occurs in the 
analysis of Figs. (3b) and (3c). Similar cancellations also appear in other regions 
of the parameter space. In Fig. (3d) the effect of cancellations in dn is shown 
as a function of mo for three sets of input data for the case when 6'^ and a^o 
have opposite signs. In each case there are large cancellations which lead to the 
appearance of minima. Aside from the reduction of the EDMs by cancellations, 
there are regions of the parameter space where kinematical suppressions occur. 
An example of this is the reduction of the EDMs when tan/3 becomes small as 
may be seen in Fig.(2f). A kinematical suppression of the EDMs can also occur 
if mi/mQ << 1. As one can see from Eqs.(19) and (29) that in this case the 
quark and the lepton EDMs are kinematically suppressed. A suppression of this 
type appears to arise in supersymmetric models with anomalous U(l) mediated 
supersymmetry breaking [25| . 



Finally, we exhibit in Fig.(3e) the excluded regions in the mo-mi/2 plane under 
the constraints given by the current experimental upper limits on and dn- The 
regions between the axes and the curves are the excluded regions in Fig.(3e). The 
analysis of Fig. (3e) shows the dramatic effect of the destructive interference on the 
allowed and the disallowed region in the mass plot. One finds that destructive in- 
terference softens significantly the stringent constraints on mo and mi/2- Thus the 
excluded region in the mo — mi/2 plane for the destructive interference case is much 
smaller than for the constructive interference case. The analysis of Fig.(3e) illus- 
trates another interesting phenomenon alluded to earlier. One finds from Fig.(3e) 
that for the constructive interference case the dn experimental constraint is the 
more severe one as it eliminates a larger part of the parameter space, while for the 
destructive interference case the d^ experimental constraint is the more severe one 
as it excludes a larger part of the parameter space in the mo-mi/2 plane in this 
case. 

In the above we have discussed cancellations which can result in a drastic 
reduction for the case of the neutron edm. There can also be cancellations for the 
case of the electron edm between the chargino and the neutralino contributions. 
For comparable sizes of 6^ and oaq, the chargino contribution is much larger than 
the neutralino contribution and cancellation is not very effective. However, more 
significant cancellations can occur for very small values of 6'^ and for moderate 
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values of since in this case the contribution from the chargino exchange and 
the neutrahno exchange become comparable. 

6 Conclusion 

In this paper we have presented an analysis of the EDM of the neutron and of 
the charged leptons within the framework of supergravity grand unification un- 
der the constraint of radiative breaking of the electroweak symmetry. All the 
supersymmetric one-loop contributions to the EDMs were analyzed taking care 
of their relative signs. For the neutron we considered also the contributions from 
the chromoelectric and from the purely gluonic operators. One finds that there 
exist significant regions of the parameter space where cancellations occur among 
the different contributions for the case of the neutron electric dipole moment. In 
these regions the neutron EDM undergoes a significant reduction and the current 
experimental limits are consistent in these regions with CP violating phases which 
are not too small and with a SUSY mass spectrum which satisfies the naturalness 
constraint. One also finds that regions of the parameter space exist where the de- 
structive interference between the different components can reduce the magnitude 
of the neutron EDM even below the magnitude of the electron EDM. 

The nature of interference, i.e., constructive vs destructive, for the neutron 
EDM determines which of the two experimental upper limit constraints, i.e., the 
upper limit on the neutron EDM, or the upper limit on the electron EDM, will 
constitute the more stringent constraint. For the case of constructive interference 
for dn, it is the experimental upper limit on d„ itself which is found to be generally 
more stringent constraint than the upper limit constraint on de- However, for the 
destructive interference case for dn, one finds that it is generally the upper limit 
constraint on de which becomes the more stringent constraint. 

As mentioned already the previously known mechanisms for the suppression of 
the neutron EDM in SUSY theories consist of suppression either by a fine tuning 
using small phases or by a choice of a heavy SUSY spectrum. We have pointed out 
a third possibility, i.e., that of internal cancellations, which naturally suppress the 
neutron EDM without the necessity of either having very small phases or having 
an excessively heavy SUSY spectrum. The cancellations that occur do not consti- 
tute a fine tuning. Rather, one finds that such cancellations occur naturally over 
a large part of the parameter space, and in some regions the cancellations become 
exceptionally large. This result has important implications for the discovery of 
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super symmetric particles. With the cancellation mechanism the SUSY spectrum 
within the current naturalness limits can be consistent with the present EDM 
experimental constraints without the finctuning of phases, and such a spectrum 
should still be within reach of the LHC. At the same time one also expects that 
if SUSY phases are indeed 0(1 — 10^^) and the SUSY spectrum lies in the usual 
naturalness limit of 0(1) TeV, then with the suppression of the neutron EMD via 
the cancellation mechanism the neutron and the electron EDMs should become 
visible with improvements of 0(10) in the sensitivity of the EDM experiments. 
Finally wc point out that although our analysis has been done in the framework of 
supcrgravity unification with soft SUSY breaking sector parametrized by six pa- 
rameters (including two CP violating phases), the mechanism of internal cancella- 
tions pointed out in this paper which can suppresss the cdms should be applicable 
to a wider class of models such as models with non-universal soft SUSY breaking. 
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8 Appendix A 

The squark(mass)^ matrix 

^^-U|i M|J' ^^^^ 
is hermitian and can be diagonalized by the unitary transformation 

DlMlD, = diag(M|i , Mj) (48) 

where one parametrizes Dg so that 

_ / cos^ -sin^e-^'^A 
^"-[sm'-^e^^^ cos'-i J' ^^^^ 

Here M?2i = \M?2i\e^'^'' and we choose the range of 9q so that ^ < Oq < ^ where 
tan^g = ^ . The eigenvalues M--^ and M~2 can be determined directly from 
Eq.(48) or from the roots 

= l(Ml, + M|22)(+)(-)^[(M|,, - + MMlilf^. (50) 
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The (+) in Eq.(50) corresponds to choosing the structure of the matrix M? so that 
for M?ii > M|2 one has M?^ > M| and vice versa. For our choice of the 9q range 
one has 

_ 2mq\Aqmo-ix*Rq\ 

where it!„ = cot P and — tan /3. Further 

g.^, ^ mo|Ag|sinag + |^|sin^^i?g 

' \moAg- ix*Rg\ ■ ^ ^ 

Using the above we get 

Im(r;i) = -Im(rf ) = lsin0,sin^, (53) 



where 





lAqUlo - H*Rq\ 









(54) 



the [+(-)] in Eq.(54) depends on whether M^^^-Ml^^ is [> 0(< 0)]. Thus Eq.(53) 
gives 

Im(rJ^) = J^r^^JfT('>T^o\Ag\ sinag + sin^^i?^), (55) 

which holds quite generally, i.e., for the case where > M~2 and for the case 
where < M~2- Thus the gluino contribution to the EDM of the quark is given 

by 

<We = ^'"^0^-(r;')[4B(^) - ^B(^)]. (56) 

One may expand the right hand side of Eq.(56) around the average squark mass. 
Defining M|=(M? +M|J/2, and expanding in the difference (M?-M| ),one finds 
in the lowest approximation 



df-giuino/(^ - ^^9<9?^(^o|^g| siu + siu ^^i?,) 



2 2 2 

mi , mi , mi 



where D(r) is given by 

1 , , Grlnr 



D{r) = ^^^^(5 + r + 2/nr + j^). (58) 



16 



As mentioned already currently there is some confusion in the literature re- 
garding the sign of the gluino contribution to the electric dipole operator p, 0]. We 
first compare our results with those of ref. The analysis of 0] corresponds to 
neglecting the D term in Eq.(57) and using dn — ^dd which gives 



2 



dn — Sa^ |.mo|Ad| sin + sin 6^^ tan/5^^^ wi^ 
e 



^ — (59) 

a a 



This result then agrees both in sign and in magnitude with Eq.(3) of ref. 0]. To 
compare with the result of ref. |0 we switch the sign of the mu term in their Eq.(6) 
(see e.g., ref. [^) and find that our Eq.(19) differs from Eq.(14) of ref. |^ by an 
overall minus sign. A comparison of the chargino and the neutralino contribu- 
tions with those of ref. |Q is more involved since the chargino (and the neutralino) 
mass matrices are different in the two works. This difference arises because after 
SU{2)l X U{\)y breaking to U{\)em^ the authors of ref.p expand the potential 
around the VEV so that Hi ^ Hi — {Hi) instead of Hi ^ Hi + {Hi) and they use in 
the chargino case Mq instead of Mc as is conventionally done[^ . Thus to compare 
with their expressions we have to do the transformation: Vij — > Cjiji, f/,*- C^ji, 
D S and X N. After that, and assuming the conventional expansion around 
the VEV, we go to their convention by the transformation Hf —>■ —Hf to find that 
we have the same overall sign in the case of the chargino exchange but the sign 
of the Kf/ term in the brackets in their Eq.(lO) should be positive. In the case 
of neutralino exchange our result differs from their Eq.(12) by an overall sign and 
further we find that the second term in the last bracket of their Eq.(13) (the term 
which begins with -Kf) should have an opposite sign. 



9 Appendix B 

The chargino matrix Mc is not hermitian, not symmetric and not real because /i 
is complex. Mq is diagonalized using the biunitary transformation 

U'*McV-^ = Md (60) 

where U' and V are hermitian and Mjy is a diagonal matrix but not yet real. U' 
and V satisfy the relation 

V{mIMc)V-^ = diag(|m^+p, |m^+|2) = U'*{McMI){U'*)-' (61) 
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We may parametrize U' so that 

cosf sm-e-^N , . 

^ - ^_sin|e-^i cosf j' ^^^^ 

where 

. 2\/^mvv^[m2 cos^ /3 + IwP sin^ + |u|m2 sin2/3cos^j5 

tan Oi — j — ^ (63) 

"^2 — 1/^1 ~ 2m^ cos 2/3 

and 

. , |//|sing^sin/3 

tan 01 = — p-7 — : — - (64) 

m2 cos p + |//| cos sm p 

Similarly we parametrize V so that 

/ cosf sinfe-^-^n , . 

^=Usi„|>^ cosf j' 

where 

^ 2v^mvK[^7i2sin^/5 + l/iPcos^/5 + |/i|m2sin2/3cos^j5 

= mi-|/.p + 2m|,cos2/3 ^^'^ 

and 

tan . -Hsin^.cos^ ^^^^ 



We wish to choose the phases of U' and V so that the elements of Md will be 
positive. Thus we define U — H x U' where 

such that 

where 71 and 72 are the phases of the diagonal elements in Eq.(60). Our choice of 
the signs and the roots is such that 

1, 



1 



[(m^ - l/xp)^ + 4m^ cos^ 2/5 + 4m^ 

(m2 + |//|^ + 2m2|//|cos6'^sin2/?)]^ (70) 

where the sign chosen is such that m + < ih^+ if 

m2 < Ia*!^ + 2m^ cos 2/?. (71) 

For the neutralino matrix, the eigenvalues and the diagonalizing matrix X must 
be estimated numerically. 
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Figure Captions 



Fig. (la): One loop diagram contributing to the electric dilpole operator where the 
external photon line ends on an exchanged chargino line labelled by xf in the loop. 

Fig. (lb): One loop diagram contributing to the electric dipole operator where the 
external photon line ends on an exchanged squark (slepton) line represented by 
qk{lk) on the internal line. 

Fig. (2a): Plot of the magnitude of the electron EDM as a function of mo when 
\Aq\ — 1.0, tan (3 — 3.0 and uao — O^o — ^ for different values of mi/2. The dotted 
curve is for mi/2 = 500 GeV, the solid curve for mi/2 = 600 GeV, and the dashed 
curve is for mi/2 = 700 GeV. 

Fig. (2b): Plot of the magnitude of the neutron EDM as a function of mo for the 
same parameters as in Fig. 2(a). 

Fig. (2c): Plot of the magnitude of the electron EDM as a function of mi/2 when 
\Aq\ — 1.0, tan/3 = 3.0 and a^o — — ^ for different values of mo. The dotted 
curve is for mo = 500 GeV, the solid curve for mo = 1000 GeV and the dashed 
curve is for mo = 1500 GeV. 

Fig. (2d): Plot of the magnitude of the neutron EDM as a function of mi/2 for the 
same parameters as in Fig. 2(c). 

Fig. (2e): Plot of the magnitudes of the neutron, the electron and the muon EDMs 
as a function of 9^^ for the case when \Ao\ — 1.0, tan /3 = 3.0, ckao = "^o = 1000 
GeV and mi/2 = 500 GeV. 

Fig. (2f): Plot of the magnitudes of the neutron, the electron and the muon EDMs 
as a function of tan/5 for the case when \Ao\ — 1.0, aAo — d^o — ^, mo = 2000 
GeV and mi/2 = 500 GeV. 

Fig. (3a): Plot of the magnitudes of the electric dipole contribution, of the color 
dipole contribution, of the purely gluonic contribution, and of the total neutron 
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EDM as a function of \Aq\ for the case when tan/3 = 3, O^o = aAo = |, 
m-g = 800 GeV (mi/2 = 281.2 GeV) and mo = 1500 GeV. 

Fig. (3b): Same as Fig. (3a) except that aAo = ~f • 

Fig. (3c): Plot of the magnitudes of the neutron, the electron and the muon EDMs 
as a function of mi/2 for the case when tan/3 = 3, 6^q = aAo = ±f , itt-o = 800 
Gev and \Ao\ = 2.6. The curve 1 (dotted) is for the case when aAo = f and curve 
2 (solid) is for case when aAo = ~f • 

Fig. (3d): Plot of the magnitudes of the neutron EDM as a function of niQ for three 
cases when O^^q = ^, aAo = — |, and tan/j = 3. The data for the other SUSY 
parameters is as follows: |Ao|=2.5, m^=500 GeV for curve 1, |74o|=2.0, mg=500 
GeV for curve 2, and |Ao|=2.5, mg=600 GeV for curve 3. 

Fig. (3e): The excluded regions in the mo — mi/2 plane of the minimal SUGRA 
model under the experimental constraints of Eqs.(44) and (45) when \Ao\ = 1.4, 
tan /3 = 3.0, 6'^o = ^ and aAo = ±f • The neutron EDM curves are solid with n(±) 
corresponding to aAo = =tf > and the electron EDM curve is dotted and labelled 
e(+, — ). The excluded regions of the parameter space lie between the axes and 
the curves. 
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